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Abstract 

We consider to detect the electron spin of a doped atom, i.e., a nitrogen or a phosphorus, 
caged in a fullerene by currently available technique of the scanning tunneling microscope (STM), 
which actually corresponds to the readout of a qubit in the fullerene-based quantum computing. 
Under the conditions of polarized STM current and Coulomb blockade, we investigate the tunneling 
matrix elements involving the exchange coupling between the tunneling polarized electrons and the 
encapsulated polarized electron, and calculate the variation of the tunneling current with respect to 
different orientations of the encapsulated electron spin. The experimental feasibility of our scheme 
is discussed under the consideration of some imperfect factors. 



PACS numbers: 03.67.-a, 61.48.-fc, 85.65.+h 



I. INTRODUCTION 



Quantum information processing (QIP) based on solid-state materials has attracted much 
attention over past few years due to the potential of scalability. Since the qubits are encoded 
in individual electron spins in most of the solid-state QIP schemes, how to efficiently detect a 
single electron spin has recently become a focus. Although there is no fundamental restriction 
for the single spin detection, it seems that, to do the qubit readout very well, we are expecting 
further development of detection technique with sufficient spatial and temporal resolution. 

Some preliminary experiments have reached the single spin level so far. For example, 
Scanning- Tunneling- Microscope (STM) electron spin resonance (ESR) has demonstrated 
single molecule ESR spectroscopy of iron impurities in Silicon yj, although theoretical work 
remains to clarify for the best description of the effects A very recent STM experiment 
under low temperature and with a high magnetic field has shown the spin flip spectra of a 
single manganese atom 3], and the magnetic resonance force microscope, which utilizes a 
cantilever driven by spins oscillating in resonance, has successfully detected a single surface 
electron spin j^J. Moreover, it has been reported that micro-superconducting quantum 
interference devices are capable of distinguishing large spin difference (Arris ~ 30) Jsj]. 

There have been some theoretical proposals for single spin detection in the candidate 
systems of QIP jslQ]. The present work will focus on the detection of a single electron spin 
of an encapsulated atom inside a fullerene (i.e., C 60 ), which corresponds to the readout of a 
qubit inperformance of quantum computing based on the endohedral fullerenes N@Cqq or 
P@Cq 0, E| • The electron spin inside the fullerene plays the role of a qubit or an auxiliary 
qubit. As the electrons of the doped atoms cannot escape the C 60 cage while preserving 
their spin states, we must develope special proposals for the spin detection, different from 
in 0, Q|. In fact, there have been some schemes in this respect, for example, by using a 
modified single molecule transistor (SMT) Q], by a nanomagnet molecule Fe§ 11 1, and by 
a shuttling device However, no relevant experiment has been achieved yet. 

In the present paper, we investigate the possibility to detect a single electron spin inside 
the fullerene by STM. STM is a mature technique available to manipulate single atoms with 
high precision at low temperature . Although no reliable evidence has been shown to 
achieve a sin£ 
proposals [if 



e spin detection, there have been some experiments P, ^ and theoretical 
12, ^ regarding single spin problems with the STM. Our scheme is different 
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from those proposals in following points: First, the electron spin to be detected, in our 
case, is assumed to be well polarized along an always-on magnetic field, while we don't 
know whether the spin is up- or down- polarized. This corresponds to the final result of 
a real quantum computing, i.e., the qubits returning to product states from entanglement 
and superposition after the elaborately designed logic gates have been carried out. So our 
detection is to know the orientation of the spin polarization. Second, as the tunneling 
through the air is different from that through a fullerene, our treatment is to separate the 
whole tunneling process into three consecutive steps. Third, due to the C 60 cage, the local 
electron spin (i.e., the spin to be detected inside the cage) is well protected. As a result, 



different from the case in 
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181 ] with the electron spin on the surface, the main error 



in our case is from the vibration of the fullerene, instead of the spin scattering, due to the 
tunneling electrons. 

The conditions we employ in the treatment below include Coulomb blockade 119] and 

n 

the polarized current in STM [20], both of which have been achieved experimentally. The 
Coulomb blockade restricts the fullerene to be charged by no more than one electron, which 
means an electron from the STM tip could jump on the Ceo only after the previous electron 
sitting on the Cqo has jumped away. This has been demonstrated in a recent experiment 
jl^ for SMT. So we will suppose that this Coulomb blockade works throughout our scheme 
under a suitable bias voltage of the STM. The polarized current in STM implies that the 
electrons going out of the STM tip are well polarized, i.e., up or down-polarized. It is 
reported in 20] that the high-quality polarized current is already available in STM. So we 
suppose below that the tunneling electrons out of the STM tip are up-polarized. In the 
presence of a magnetic field, the tunneling electrons will couple to the caged electron by 
exchange interaction due to the spin degrees of freedom involved 2lj]. Moreover, the total 
electron spin of the doped atom N or P are S = 3/2 with four Zeeman levels | ±3/2) and 
| ±1/2) It has been demonstrated that quantum gating can be performed independently 



with electron spin states | ± 3/2) or | ± 1/2) [22|. For simplicity, however, we will focus 
below on the discussion about the caged electron spins | ± 1/2). The case regarding | ± 3/2) 
could be obtained simply by enlarging the variation of the tunneling current by three times. 

In the next section, we will explore the spin-dependent tunneling matrix elements. Then 
a specific calculation of the tunneling will be made in Sec III which also includes discussions 
about the experimental imperfection. The last section is for the conclusion. 



3 



II. THE SPIN-DEPENDENT TUNNELING 



To clarify our description below, we simplify the system shown in Fig. 1 to be a doped 
fullerene sandwiched by two leads L and R, as shown in Fig. 2. So the Hamiltonian is, 

H = Y1 e \kc{ Al C XA l + H c + ^A,k,l a n,l + ^A c A,fc,2 a n,j)> (!) 

\,k,l X,k,n,l 

where c\ k l (c\.k t i) creates (annihilates) an electron in the lead X = L,R, with k being the 
momentum of the electron, I = ±1 for (up/down) spin polarization, and are energies 
regarding Fermi energy. The operator a nl {a n ,i) is related to the electron on the fullerene 
in the orbital n. t\ is the tunneling coefficient regarding A = L or R. H c accounts for the 
fullerene including the energy of itself and the Coulomb blockade term %ak)(^2k a \ a k ~ 
1) as well as the exchange coupling — JS ■ a with a and S the spin degrees of freedom of the 
electrons outside and inside the fullerene, respectively. Since our interest is in the tunne ling 
modification due to the electrons' exchange interaction, instead of the electro-phonon 23[ 
or spin-phonon resonances, we may follow the idea in jl^l . [nl Q] and focus on the spin 
dependent tunneling matrix elements. To this end, we separate the tunneling into three 
steps (See Fig. 2): The first step is from the left lead to the fullerene, the second one is for 
tunneling through the fullerene, and the third accounts for the tunneling from the fullerene 
to the right lead. 

Both the first and the third steps can be modeled as electronic tunneling through ca- 
pacitors with $j = Widi where i = 1,3, $j is the work function, dj is the tunneling gap, 
and Wi represents the constant regarding the free space permittivity, the dielectric constant 
and the area of the lead. By using WKB approximation, we could reach the spin-dependent 
tunneling matrix elements as below [24 1, 



and 



T l0 cexp{-W^/ V $ i- JS-adx}, (2) 



I <• >l \/^J V $3 -JS-adx}, (3) 

where J is the strength for the exchange coupling between the tunneling electrons and the 
caged electron due to the spin degrees of freedom of the electrons, which is proportional to 
the integral below 

e 2 

*t(ri)^(r 2 ) — ^ l {r 2 )m 2 {r 1 )dr 1 dr 2 (4) 

T 12 




with ^i(ri) and ty 2 (rj) the wavefunctions of the electrons tunneling and caged at the position 
Ti and rj, respectively, and r\ 2 — \t\ — t 2 \. For simplicity, we have omitted in Eqs. (2) and 
(3) the unimportant constants regarding the spin independent tunneling rate in the absence 
of J. Please note, we have also made an assumption here that, due to the Coulomb blockade 
and the low temperature limits, each electron from the left lead would jump to the ground 
state of Cq and the tunneling amplitude is mainly dependent on the exchange interaction, 



instead of the electronic momenta. Since both spins of the tunneling anc 
are well polarized, we may simplify JS ■ a to be Ja z with J being 1 meV 
straightforward algebra for Eq. (2) yields corresponding tunneling rates, 



caged electrons 



2l|. As a result, 



2 /8mW lu , J 3/2 _ _J_^ 3/2] 



eXp| 3 V h 2 1 ( 2WW 1 ' 
where the signs '±' correspond to the up/down polarization of the caged electron spin. 
We have approximately omitted in above equation the high-order terms regarding J/Wi 
because they are comparatively negligible (i.e., J/W\ ~ xi/1000 to the numbers shown 
later). Similarly, we have 



2 8mW 3 3/2f 3 J 3/2 , 3 J 

2 /8m^ 3 3/2. , 3J 2 l8mW 3 3/2 , 37 



W eXp{ -3V^ X3 (1± 2H^ )}eXp{ 3V^ X2 (1± 2W^ )} - 

The second step can be considered as a tunneling through a constant potential barrier. 
This is because the fullerene could be modeled as a spherical capacitor. Due to the equipo- 
tential surface of the spherical capacitor, the work function $ 2 should be a constant, which 
yields the tunneling rate as 



T 2 ocex P {-y'^($ 2 ±J)}. (5) 
Since J is smaller than $j by 1000 times, we may write the terms regarding J as 

exp{±A} = cosh(A) ± sinh(A) « 1 ± A, (6) 

where D\ = J \/8mx\ /(Wih 2 ) is from Ti, D 2 = Jy/8 mx 3 /(W 3 h 2 ) and D 3 = 



J^8mx 2 / (W 3 h 2 ) are from T 3 , and D 4 = \ / 2mJd 2 /y / <& 2 h 2 is from T 2 . 



Since the tunneling current is proportional to the multiplication of the tunneling rates 
above, we may simply write the tunneling current as 
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J = / J](1±A), (7) 

i=i 

where 

i.e^-^'Vpf-^l-f-fll^l-^A (8) 

is the spin independent tunneling current, which could be detected in the absence of the 
magnetic field. In the following calculation, we assume Iq to be 1 nA if d\ = 0.23 nm. 



III. DISCUSSION 



From [19| we know that, d 2 = 0.7 nm is the size of the fullerene, and = 0.27 nm is 
due to van der Waals interaction between the fullerene and the substrate. We suppose that 
is 2 eV between x = and x — x\, $2 = 1 eV is constant, and $3(2) = 1 eV is from 
x = X2 to X3. Then we can obtain by direct calculation that AI± = I± — Iq = ±16 pA, 
where I± corresponds, respectively, to ± in Eq. (7). The current of 16 pA is detectable 
with the present STM technology. Slightly changing d\, we may estimate AI = AI + — A/_, 
which changes drastically with respect to di, as shown in Fig. 3. So it is evident that the 
difference between spin-up and spin-down of the caged electron is distinguishable from the 
tunneling current as long as d\ is around 0.24 nm. 

As mentioned above, because the local spin in our case is pretected by the fullerene cage, 
the main source of error is the vibration of the fullerene due to the tunneling electrons, 
which yields a 8 h lft „ = 3 p m e,— to the ene rgy variation 5 m eV Q. Became this 
vibrational degree of freedom is far detuned from other characteristic frequencies, and also 
because there is no evidence of vibration-spin coupling, we did not consider this vibration in 
above treatment. The tunneling current, however, is very sensitive to the distance variation. 
So we have to strictly assess the influence due to the position shift of the fullerene on the 
tunneling current. To this end, we consider the tunneling gaps regarding the first and the 
third steps to be changed time-dependently, while the second step remains unchanged. So 
we have the modification of the tunneling rates, 



2 SmW, 3/2 3A ± 3J 
Tl K eXp{ ~3 V (1 + 2^ ± 2W } ' (9) 



T 3 oc M^&p ± -i^) } exp { ?.S I f(l - ^ ± «_)>. (10) 
3 P1 3 V ft 2 3 2t*W _l3V /I 2 2 V 2x 2 2W 3 x 2 Jf V ; 

Suppose A = (8/2) cos(u;i) with a; = 1 THz |l9j . Like in jl6|, we introduce the average 
current over the time T, i.e., (I) = (1/iV) J^Zi I ( t i)^ where the 

sum from i = 1 to N is 

over the number of the tunneling electrons in the time T, and is proportional to the 
tunneling rate for each electron. So (I) can be written as 



(I) = IoT\(l±D i )[l±D 1 + (5/2x 1 )J—-± / (l/T)cos(ujt)dt][l±D 3 (11) 

i=2A V 11 J 



8mW 3 



T 



-(6/2x 2 )y—jp-± J (l/T)cos(ut)dt}. 

As the terms regarding 5/x\(2) are comparable to -Di(3), we have to eliminate the vibrational 
influence. To this end, we may detect the current during a time period T = kir/u with 
k — 1, 2, • • •, in order to average out the terms regarding 5. This would in principle make 
the fullerene vibration, due to the tunneling of the electrons, of negligible effect on the 
original current in our observation. However, due to the intrinstically stocastic property of 
the electron tunneling, the efficiency of above trick depends on the statistical deviation a 
of a normal distribution (1/(a/27t<j) exp [— (t — to) 2 / (2a) 2 ], where to is the mean tunneling 
time of each electron we desire Q]. To carry out our scheme to the best, we require a to 
be as small as possible. 

Besides, there are other possible imperfection in implementing our scheme. In a real 
quantum computing, the qubits are not usually well polarized at the end of the computing 
operations. For example, the Cgo-caged electron-spin would probably be in superposition 
F\ |) + G\ [) (F > G or F < G) due to some unpredictable error sources in quantum 
computing. The estimated current dispersion due to this imperfect polarization is e.g., 
F> G, 



(A/ 2 ) 2G * 

"7I?>' (12) 



i=l 

in which G is much smaller than 1 and iV is of the order of hundred in an observation during 
a period of T =100 ns. So this current dispersion is too small to affect our spin detection. 
In addition, the feedback effect of the tunneling electrons on the caged spin should also 
be addressed. This estimate could be done by Fermi golden rule from the second-order 



perturbation 



12 1, which yields the spin decay rate 1/t s = \ Ylt=i Di\ 2 /r e with l/r e the 
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tunneling rate of the electrons. As we suppose Iq = 1 nA, implying l/r e = 10 10 Hz, we 
have r s > 10~ 6 sec, which is much longer than our detection time. So this spin decay is also 
negligible. 

Normally, the tunneling current in usual use of STM is from the order of pA to nA, which 
implies that the tunneling rate of the mobile electrons is from one electron every 100 ns to 
one electron every 100 ps. We have noticed that the observation of individual electron-j ump 
is within the reach of the current technique by I-V characteristic and by dl/dV plot [25], 
and the minimum current the STM is able to distinguish is of the order of 0.1 pA 1261. 

n 

Recent experiment [27] has shown the availability of a controllable manipulation of a 
single Cgo by STM. Moreover, in the absence of the tunneling of the electrons outside the 
fullerene, the decoherence time of the caged electron-spin is about 1 second at 7° K much 
longer than our implementation time. So decoherence is neglected in our treatment above. 
While there is no experimental data for dephasing of the tunneling electrons. A recent 
work j3] has shown that the spin dephasing time of the bulk two-dimensional electron gas 
at low temperature could be 150 ns. In our case, we require the dephasing time of the 
tunneling electrons to be longer than 100 ns. Considering the data in 28j and the difference 
between the bound and unbound electrons, our requirement should be satisfied at very low 
temperature. 

Compared with a previous work using a modified SMT, the present scheme is, to some 
extent, similar, but much simpler. In spin flips of the mobile electrons by microwave 
pulses are necessary. Since the electron tunneling is intrinsically stocastic, the decoherence 
due to imperfect spin flip yields the main infidelity. In contrast, no spin flip is needed in 
the present scheme. So the dephasing of the mobile electron is not a serious problem in our 
scheme. Furthermore, a magnetic field gradient dB/dz = 4 x 10 6 T/m is essential to [10]. 
While such a magnetic field being highly stable in time and very large and homogeneous 
in spatial gradient is still challenging with the present technique. In contrast, the constant 
magnetic field required in the present scheme is fully within the reach of the current tech- 
nique. Shortly speaking, the present scheme is better than in Q], and the advantages are 
from a combination of STM and SMT: STM makes sure that the tiny variation of the cur- 
rent can be distinguished, and the Coulomb blockade in SMT experiment guarantees that 
the change of the tunneling current corresponds to the fullerene charged by only a single 
polarized electron. 



S 



IV. CONCLUSION AND ACKNOWLEDGEMENT 

In summary, we have explored the possibility to detect the electron spin of the doped 
atom in a fullerene cage by STM, assisted by the Coulomb blockade and polarized current. 
Different from previous studies, we separated the whole tunneling process into three parts. 
After presenting the spin dependent tunneling matrix elements, we have shown by specific 
calculation the feasibility of this detection with current STM technique. The detrimental 
influences from some imperfect factors have also been considered in our treatment. Since the 
STM is widely employed in atomic and molecular control in various systems, we argue that 
the main idea of the present work, i.e., the use of the STM device along with the Coulomb 
blockade and the polarized current could be in principle applied to other candidate systems 
of QIP for detection of a single electron spin. 
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The captions of the figures 

Fig. 1 Schematic of the detection of a single caged electron spin by STM, where the black 
dots mean the electrons, and a and S correspond to the spin degrees of freedom regarding the 
electrons tunneling and caged, respectively. The caged electron spin can be ±1/2 or ±3/2, 
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which pair of spins employed in quantum computing has been known before the detection. 
So what we want to detect is whether the spin is up (i.e., 1/2 or 3/2) or down (i.e., -1/2 or 
-3/^2). The gap between the fullerene and the substrate is due to van der Waals interaction 

Fig. 2 (a) Schematic of our treatment to model the spin-dependent tunneling by three 
steps, and (b) the work functions in different steps are assumed. 
Fig. 3 The current difference AI = AI + — Ai_ with respect to d\. 
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